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Abstract 

We study a class of statistical systems which simulate 3D gonihedric sys- 
tem on euclidean lattice. We have found the exact partition function of the 
3D-model and the corresponding critical indices analysing the transfer ma- 
trix K(Pi,Pf) which describes the propagation of loops on a lattice. The 
connection between 3D gonihedric system and 2D-Ising model is clearly seen. 
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1 Introduction 



In the articles the authors formulated a model of random surfaces with an action 
which is proportional to the linear size of the surface. The model has a number of 
properties which make it very close to the Feynman path integral for a point-like 
relativistic particle. In the limit when the surface degenerates into a single world 
line, the action becomes proportional to the length of the path and the classical 
equation of motion for the gonihedric string is reduced to the classical equation of 
motion for a free relativistic particle.^] 

In addition to the formulation of the theory in the continuum space the system 
allows an equivalent representation on Euclidean lattices where a surface is associ- 
ated with a collection of plaquettes @, U . In these lattice spin systems the interface 
energy coinsides with the linear-gonihedric action for random surfaces. This gives an 
opportunity for analytical investigations || [7], || 0, [14]] and numerical simulations 



[ jig , [lq , [17], [18|, [IjJ of the corresponding statistical systems. 

Additional understanding of the physical behaviour of the system comes from 
the analysis of the transfer matrix which describes the propagation of the closed 
loops-strings in time direction with an amplitude which is proportional to the sum of 
the length of the string and of the total curvature. In this article we shall study the 
physical picture of string propagation which was suggested in the transfer matrix 
approach ||. 

The partition function of the system is defined as [|TJ 

ZgonihedrUP) = E tX V{ ~0A(M) } , A(M) = £ Ay ■ |tT - Ofy], (1) 

{M} <ij> 

where Ay = a is the length of the edge < ij > which is equal to the lattice spacing a 
on the cubic lattice, is the dihedral angle between two neighbouring plaquettes 
of the singular surface M sharing a common edge < ij >, cty = 0,7c/2,n. ^ In 
([]]) {M} denote the set of closed singular surfaces on the three-dimensional toroidal 
lattice T 3 of size N x N x N. A singular surface M is a collection of plaquettes in 
T 3 such that every link is contained in 0, 2 or 4 plaquettes and every plaquette of 
the lattice T 3 can be occupied only once f|, ||. The surfaces are closed because only 
an even number of plaquettes meet at a given lattice link. The singular surfaces of 
interface which describe the states of arbitrary three-dimensional spin system can 
be viewed as a set of surfaces {M} J|, [HJ. 

In it has been proven that the partition function ([!]) can be represented in 
the form 

Z(fl)= K p (P 1 ,P 2 )---Kp(P N ,P 1 )=trK», (2) 

{P!,P 2 ,...,P N } 

where Kp(Pi, P 2 ) is the transfer matrix of size 7x7, defined as (see formulas (15), (16) 
of i) 

Kf mhedrtc (P u P 2 ) = exp{-P [fc(Pi) + 21{P X A P 2 ) + k(P 2 )} }, (3) 



lr The problems of spiky instability and the convergence of the partition function have been 
studied in § g § 

2 Usually we take a = 1. 



where Pi and P 2 are closed polygons on a two-dimensional toroidal lattice T 2 of 
size N x N and 7 is the total number of polygon- loops on a toroidal lattice T 2 . 
Closed polygons {P} = II are associated with the collection of links on T 2 with the 
restriction that only an even number of links can intersect at a given vertex of the 
lattice and that the links can be occupied only once. Q 

The transfer martix can be viewed as describing the propagation of the 
polygon-loop Pi at time r to another polygon-loop P2 at the time r + 1. The 
functional k(P) is the total curvature of the polygon-loop P which is equal to the 
number of corners of the polygon ( the vertices with self-intersection are not counted) 
and l(P) is the length of P which is equal to the number of its links. The length 
functional /(Pi A P 2 ) is defined as (see formula (12) of || ) 

/(Pi A P 2 ) = /(Pi) + Z(P 2 ) - 2/(Pi n P 2 ), (4) 

where the polygon-loop P x A P 2 = P x U P 2 \ Pi fl P 2 is a union of links Pi U P 2 
without common links Pi fl P 2 . The operation A maps two polygon- loops Pi and 
P 2 into a polygon- loop P = Pi A P 2 . Note that the operations U and fl do not 
have this property. These operations acting on a polygon-loops can produce link 
configurations which do not belong to ILR The length functional /(Pi A P 2 ) defines 
a distance between two polygon-loops Pi and P 2 . It is natural that the transition 
amplitude Kp(P\, P 2 ) (|3|) depends only on the functional which mesures the distance 
between the two subsequent polygon-loop configurations, because transition ampli- 
tude decreases when the distance (f|) between two configurations increases. This 
can be seen also from the inequality 

/(PiAP 2 ) > \l{P x )-l{P 2 )l (5) 

therefore 

K 9 ; mhed " c (P u P 2 ) < exp{-2(3 |/(P) - /(P 2 )| }. (6) 

Algebraically one can construct many functionals of that kind, but what is important 
here is that this distance functional appears naturally from the geometrical action 
([!]) of the original theory. The expression for the transition amplitude Kp(Pi, P 2 ) in 
the polygon-loop space II is very close in its form with the transfer matrix for the 
random walks 

K p (X,Y) = exp{-P\X-Y\}, (7) 

which depends only on the distance between initial and final position of the point 
patricle. Statistical mechanics of paths with curvature-dependent action is also well 
known fllj. The aim of this work is to study spectral properties of the transfer 
matrix K /3 (P 1 ,P 2 ) and critical behaviour of the statistical system (|2]). 



2 Free energy and correlation functions 

The eigenvalues of the transfer matrix Kp(P\, P 2 ) define all statistical properties of 
the system and can be found as a solution of the following integral equation in the 

3 The polygon-loops P appear as the intersection of the singular surfaces M with the planes 
between coordinate planes in T 3 || . 

4 The symmetric difference of sets P\ A P2 is an important concept in functional analysis Ca|. 



loop space LTf] 

J2 K P (P 1} P 2 ) *(P 2 ) = A(/3) *(P0, (8) 

where ^(P) is a function on loop space. The Hilbert space of compex functions 
*(P) on n will be denoted as H = L 2 (U). 

The eigenvalues define the partition function (^j) 

ZC9)=Af + ... + <, (9) 
and in the thermodynamical limit the free energy is equal to 

-PfW-g^^lnZVi). (10) 
The correlation lengths are defined by the ratios of eigenvalues Aj(/3)/A (/3) 

tU Ao(/3) 

and grow if the eigenvalues Aj(/3) approach the eigenvalue Aq(/3) at some critical 
temperature (3 C - By the Frobenius-Perron theorem Aq(/3) is simple and we have 

A (/3) > Ax(/3) > A 2 (/3) > ... (12) 

Finite time propagation amplitude of an initial loop Pj to a final loop P/ for the 
time interval t = Mj (3 can be defined as 

K(P h P f )=Av M K /3 (P i ,P 1 )...K f3 (P M _ 1 ,P f ), (13) 

{P 1 ,P 2 ,...,P M -i} 

where we have introduced natural normalization to the biggest eigenvalue Ao and 
M < N . The trace of the operator K(Pi, P/) is equal to 

TrK = (1 + (A!/A ) M + • • • + (A 7 /A ) A/ ) (14) 

and depends on the ratio Aj/A . 



3 Connection between 3D gonihedric system and 
2D-Ising model 

We consider below a transfer matrix which is less complicated than the original 
matrix @ and depends only on the distance functional /(Pi A P 2 ) 

Kp{Px, P 2 ) = exp{-2(3 ( /(P A P 2 ) - N 2 ) }. (15) 



5 We shall use the word "loop" for the " polygon- loop" . 

6 It is convinient to substruct the vacuum energy 2N 2 in the exponent. 



The largest eigenvalue appearing in the equation 



£ exp{-2/3/(Px A P 2 ) + 2f3N 2 } • *(P 2 ) = A(/3) (Pi 
{^2} 



(16) 



can be found because the corresponding eigenfunction is a constant function ^!q(P) 
1 (see Appendix). Therefore 



A = J2 exp{-2pl(P x A P 2 ) + 2/3N 2 }. 



(17) 



The sum (|17|) does not depend on P x . To prove this we note that the loop P = P X AP 2 
runs over all loops in II as P 2 runs over II, i.e. the mapping P x — > P% A P 2 is one to 
one for any P 2 . This change of the variable proves that 



A = E exp{-2/3l(P) + 2(3N 2 }, 
{P} 



so An is the partition function of the 2D-Ising ferromagnet. Indeed [19, EG, 21 



Z 2D - Ising = £ exp{-2pi(P) + 2(3N 2 } = e^/O^ 
{P} 

where f{0) is the free energy of the 2D-Ising model. Therefore 



An 



J 2D-Ising 



A^ + ... + A 



N 

2 N i 



;is) 



(19) 



(20) 
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where \ are the eigenvalues of the transfer matrix of 2D-Ising model [Ej 
Thus the largest eigenvalue of the SD-system ([15) is equal to the partition function 
of the 2D-Ising ferromagnet. 

The free energy of the 2D-Ising ferromagnet in the thermodynamical limit is 
given by gg, |2T], £2| 



2tt 



d^drj 



2 Jo (2tt) ; 



Zn[ (1 + w A ) 2 - 2(w 2 - w b )(cosf + cost?) ], 



(21) 



where w = e _/3 . Therefore the free energy of the three-dimensional system which is 
defined by the transfer matrix flT5|) is given by f(/3) and coincides with the one of 
2D-Ising ferromagnet. 

From this result we can deduce that the critical temperature of the three- 
dimensional system (p~5|) is equal to the one for the 2D-Ising ferromagnet 2f3 c = 
ln(y2 — 1), that the specific heat exponent a = and from the hyperscaling law 
vd = 2 — a that v = 2/3. 

We recall || that in the alternative approximation when the intersection term 
2k{P\ fl P 2 ) is ignored in the original transfer matrix ([3]) 

K 0( 3 (Pi, P 2 ) = exp{-(3 [k(P x A P 2 ) + 2l(P x A P 2 ) - 2N 2 ] } (22) 

the free energy can also be computed f| 

1 d£dr) 



-/3/o(/3) 



2 Jo (2tt)< 



Zn[ (l + w 4 ) 2 - WV(l-o/) + 



+ 4w 4 (l - u 2 ) cos f cos r] - 2(w 2 - 2wV + w 6 )(cos£ + cost?) ], (23) 

where u 2 = w is the contribution from the curvature term k(P) and exhibits the 
same critical behaviour as the 2D-Ising ferro magnet. Thus the original gonihedric 
system (EI), (Rl) is bounded by two close statistical systems 



PHP) < -PfgonihedriciP) < ~Pf(P), (24) 



because 



K p(Px,P 2 ) < K°T hedHc (P x , P a ) < K^(P 1 ,P 2 ) (25) 

This confirms the conjecture || that 3D gonihedric system should have statistical 
properties close to the ones of 2D-Ising ferromagnet. Earlier numerical simulations 
[ PS| , [T7|, [TEj support this dimensional "reduction". It is also consistent with the 
analytical estimate of the enthropy factor of the random surfaces on a cubic lattice 

Layer-to-layer transfer matrices for three-dimensional statistical systems, whose 
elements are the product of all Boltzmann weight functions of cubes between two 
adjacent layers have been considered in the literature [E5|, ER E3|, EJ1. Using Yang- 



Baxter and Tetrahedron equations one can compute the spectrum of the transfer 



matrix in a number of interesting cases p3l G4J. In the given case the transfer matrix 



@ has geometrical interpretation which helps to compute the spectrum. 

We have to remark that if the loops {P} are not restricted to be closed the 
system is essentially simplified. Indeed in the model where all subsets of links from 
T 2 are allowed as configurations of the system and the transfer matrix is defined by 
the same formula the model becomes trivial, the transfer matrix is the 2iV 2 th tensor 
product of lD-Ising model transfer matrix 

l-W ^\ ( 26 ) 
Using this observation one can derive inequality 

(1+e - 2/3) 2^-l(l_ e -2/3) 

Ai < . 



*->2D-Ising 

4 Loop space and eigenfunctions 

To proceed it is convenient to introduce some notation. An invariant product in II 
can be defined as 

<P X |P 2 > = Z(PiAP 2 ) - 1{P X AP 2 ) (27) 

where P = T 2 A P = T 2 \ P and T 2 is the loop which contains all links of the 
toroidal lattice T 2 ( l(T 2 ) = 2N 2 ). The invariant product ( fZTD is also equal to 



< p x \p 2 > = 2N 2 - 2 l(P x A P 2 ) 



(28) 



and is odd with respect to P-the complement operation < Pi\P 2 > = — < Pi\p2 >■ 
In particular the energy functional of the system is equal to 

E P = < | P > = 2N 2 - 2l(P) (29) 

The product ( |27|) is invariant under the simultaneous rotations of the loops Pi and 
P2, defined as 

P -> PA 5, (30) 

where 5 is an arbitrary loop 

<P 1 A5\P 2 A5> = <P 1 \P 2 >. (31) 

This group of transformations in the loop space II is Abelian because for its repe- 
sentations on H 

R 5 ^r(p) = $(PA«5), (32) 

we have 

[R Sl , R S2 ] = R] = 1. (33) 

The product < P1IP2 > is invariant also under translations. The group of transla- 
tions is defined as a rigit translation of the loop P in x and y directions on T 2 by 
same units of lattice spacing a 

P — ► P + CL x e x + a y e y . (34) 

Together these two groups form a Nonabelian group which acts on the loop space 

n. 

With this notation the transfer matrix fllBI) takes the form 



K p {P x ,P 2 ) = exp{(3 <P|P 2 > }. (35) 

and the integral equation ( jl6|) the form 

£ e P <Pi\P2> . #(p 2 ) = A (0) ^( Pl ) ( 36) 
{ft} 

Note that (06|) is invariant under rotations (|30|) . This suggests that we search for 
eigenfunctions of the operator (|35| ) in the form of invariant polynomials in H = 
L 2 (I1). These polynomials can be constructed in the form of powers of the invariant 
< P\Q > as follows: 

<(P) = 1, 
*g } (P) =<P|Q>, 
^ 2) (P)=<P|Q> 2 - P2 , 

¥"\p)=<p\q> 3 - p 3 <p|g>, 



^\P)=<P\Q> n - Pn <P\Q> n ~ 2 -... (37) 



where the coefficients 



E { P } <Q|P> 2 E { P } <Q|P> 4 

P2 " W ' P3 " E { P> <0|P> 2 ' (38) 

and so on are chosen so that Wq is orthogonal to if i ^ j. The set of functions 

{*q (P)} in H form an invariant subset H n of the level n. The index Q numerates 
functions inside the level. We introduce general Legandre loop polynomials as 

L n (x) =x n - p n x n ~ 2 - (39) 

so L n (< P\Q >) = ^/q\p). This set of functions in H is appropriate for solving 
the integral equation (RBI) and in its form is very similar with the ones which are 



used for the random paths with curvature-dependent action [11|. 

Let us begin by proving that ^Jq\p) is the eigenfunction of (j36|). Because (see 
the next section) 

E <Qi|P> n • <P\Q 2 > =/im<Qi|Q 2 > (40) 
{p} 



where 



we have 



where 



and thus 



Mnl = E < I P > n+ \ m 1 = (41) 

£ ^ <A|P 2 >. < p 2 |g >= Al(/3) < Pi |g > (42) 
{^} 

oo on i 

A i = E V" 1 = ^E ea; P{/ 5 <0|P> } < | P > (43) 

n=0 n - AIy {P} 

Al = 2^ | A °' (44) 
For the ratio Ai/A one can get 

" (Ai/Ao) = ^ ^ In A = „(/?) (45) 

which is the internal energy of the 2D-Ising system. Thus the second eigenvalue of 
the three-dimensional system coincides with internal energy of the two-dimensional 
Ising system. 

Using equation ( [40"D one can also address the question of degeneracy of the eigen- 
values. The scalar product of the loop functions can be defined as 

*Qi • *Q 2 = E ^Qi( P ) ^Q 2 ( p ) (46) 
{P} 

then for the zero and first level functions ( |3"7| ) we have 

•a ■ »a = o, 



= < QilQa > - (47) 
The rank of the matrix Gq 1 q 2 =< Qi\Q 2 > defines the number ei of linearly 

independent functions on the first level Hi and thus the degeneracy of the eigenvalue 
Ai f\ The number €\ is bigger or equal to the number of functions ^q] which are 
orthogonal to each other. The orthogonality condition follows from (HTJ) 

<Qi\Qj> = i,j = l,2,....,2iV 2 

The last equation can be rewritten also in the form 

l(Qi A Qj) = N 2 

and its solutions provide 2N 2 linearly independent functions on the first level. The 
solutions will be presented in a separate place. 

In the following we will search for approximate solutions since the sums 

E <Qi\P> n - <P\Q 2 >\ k>2 
{P} 

depend not only on invariant < Qi\Q 2 >> but also on the shape of Qi and Q 2 . 
However the sum depends only weakly on the shape of Q's and the leading behaviour 
involves only the length. In this approximation we have (see also the next section) 

< Qi\P > n ■ < P\Q 2 > 2 = /x„2 < Qi\Qa > 2 + /i«o (48) 

{P} 

where 
and 

VnO = -TZjJl E(<0|P> n+2 -(2iV 2 ) 2 <0|P>") 

[AIM ) - p 2 | P } 

so that the next eigenvalue can be found in the same way 

A 2 = (9N2 i - E ™ V { P < I P > } [< I P > 2 -p 2 ), (49) 

and for the third one we have 
A 3 = , OAm3 1 , OAr2 , E exp{/5(<0|P>)}[<0|P> 3 -p 3 <0|P>], (50) 



(2N 2 ) 3 - p 3 (2N 2 ) 



{P} 



and so on. All eigenvalues can be expressed through the generalysed Legandre 
polynomials d39|) 



An = Ln{ l N ^ E e MP <0|P> }L n (<0|P>), (51) 



{P} 



7 The degeneracy of the zero level is one cq = 1. 



or as corresponding derivative of the Ising partition function 

AM = L^m L ^ (52) 

This completes the computation of the eigenf unctions and eigenvalues of the integral 
equation (0), (0). 



5 Loop correlation functions 

The equality ( flOf) can be proven by using the fact that the l.h.s. of (]40[) is equal to 
X <0\P> n -<P\Q> = ]T <o\p > n+1 -2l(Q)J2 <0|P> n + 

{P} {P} {P} 

+ 4^ <0\P > n ■ l(Pf]Q). (53) 
{P} 

The last term in ( |53"D is a linear function of 1{Q). If Q = Q\ U Q2 and Qi H Q2 = 
then /(P n Q) = l(P n Qi) + Z(P n Q 2 ) and we have 

£ <0|P>"-/(Png) = tj n iZ(Q). (54) 



The normalization constant can be computed at the point Q = T 2 

1 

2iV2 



Vm = tt^E <0|P>"/(P). (55) 



{P} 

The linearity of the last term in ( |53~D with respect to l(Q) follows also from the 
decomposition of the functional l(P C\ Q) into the sum over the links of Q 

/(PnQ) = $>(P) 

where Xi{P) = 1 if i G P and zero otherwise. By homogenuity 

E <o\p> n - Xl (P) 
{p} 

is independent of i. Hence the summation over % e Q gives the expected result fl54|). 
To proceed we have to compute the sum 

X; <gi|p> n -<p|g 2 > 2 , (56) 



which by the transformation similar to (|53|) can be expressed in terms of the two-set 
correlation function 

<0|P> n • l(PnQ 1 )-l(Pf]Q 2 ). (57) 



The last correlation function is a sum over i e Q\ and j G Q 2 of the two-link 
correlation function 

£ <0\P> n - X i(P)Xj(P)- (58) 
{P} 

In the approximation when the two-link correlation function ([58"D factorises, then 
the correlation function Q57| ) is equal to 

Vn2 KQ1KQ2) + VnO KQl n Q 2 ) 

and we recover the expression (|48|). 

For the high powers we have to compute a many-set correlation function 

£ <o|P>"-]pOPnQ a ) = Vna l[i(Q a ) + ••• (59) 

{P} a a 

which can be expressed as a function of lengths in the approximation when many-link 
correlation function 

£ <0\P> n -UP)xAP)Xk(P)---, (60) 

factorises. Because the link variable Xi{P) is a quadratic function of Ising spins 
it follows that many-link correlation function (BO) is nothing else than many-spin 
correlation function of the 2D-Ising model. 
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7 Appendix 

It is elementary to verify that the function 1{P\ A P2) (|4]) on II x II is a distance 
functional p(P u P 2 ) = l(P 1 AP 2 ) and that if Pi, P 2 , P e IT are such that l(P 1 AP 2 ) = 
/(Pi A P) + l(P A P 2 ), then P 1 n P 2 C P C P l U P 2 . Indeed we have 

l(P) = 1{P n Pi) - 2/(Pi n P 2 ) + /(P 2 n P), 

hence /(P\Pi U P 2 ) = /(P n Pi H P 2 ) - /(Pi n P 2 ) from which we conclude that 
P C Pi U P 2 and P x H P 2 C P. 

We complete our arguments by showing that all eigenvalues of the normalized 
transfer matrix Aq 1 Kp(Pi, P 2 ) are smaller than or equal to 1. In order to prove 
this statement let \p De a "delta function" in II, i.e. Xp(Q) — 1 if P — Q an d 
Xp{Q) = otherwise. The set of functions {xp} form a basis in H = L 2 (U). Note 
that for the propagator (p~3|) 

if(Pi,P/) =< Xp, \^ M Kf\ XPf > 



we have (1131 



K(Pi,P f ) = A ~ M Yl K f3 (P i ,P 1 )K p (P 1 ,P 2 )-.-K /3 (P M _ 1 ,P f ) < 

{P 1 ,P 2 ,...,P M -i} 

< A M Yl K p {P i ,P 1 )Kp(P 2 ,P 3 )---K f) (P M „ 1 ,P f ). 

{ft,ft,...,PAr-i} 
We can sum successively over P/s to obtain 

K(P,P/) < Aq 1 

for any Pj, Pj provided M > 1. If K(Pi, Pf) had an eigenvalue A > 1 then K(Pi, Pf) 
would grow as A M as M — > oo for same P/. Since i ; T j a(Pi,P 2 ) > it follows from 
the Frobenius-Perron theorem that 1 is a simple eigenvalue. 
It follows that for M > 1, #(P, P) < A \ thus 

A (/3) 

where the mass m(/3) is defined as 

K(P, P) « exp{ - M m(/3) } 

for sufficiently large M. From the other side if in the sum which defines K(P, P) we 
restrict all intermediate loops to be P then Ag M (/3) < K(P,P) and thus 

'^1 < m(p) < to Ao(/3) . 

The operator K(P,P) is bounded from below also by positive number I/7. In- 
deed, because 

K(P,P) =< XP \A M K^\ XP >=Y< Xp |Ao M ^ M | * W >< * (n) \Xp >= 

n=0 

= E(^) M l<xp|* (n) >l 2 > I < XP I * (0) > I 2 . 

n=0 A 

Since for the normalized constant function on II 

< X p I * (0) >= VV7, 
where 7 is total number of loops on a toroidal lattice T 2 , thus 

l/ 7 < K(P,P). 

For the mass m(/3) one can get 
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